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Abstract 


Solar radiation that reaches the Earth's surface has a significant impact on human health, 
agricultural productivity, terrestrial, aquatic ecosystems, and air quality. The difficulties from 
shortage of experimental field data motivate the researchers to rely on other database such as 
satellites and reanalysis radiation data. In the present work, the integrated daily solar ultraviolet- 
visible radiation (total integral of ultraviolet-visible radiation from 200 to 440 nm) was extracted 
from the ERAS reanalysis database during the historical period of 1979-2020. Moreover, Total 
Column Ozone (TCO), and the Percentage of Cloudiness (CF) were obtained from the same 
database. The study area was divided into eight clusters by Analytic Hierarchy Process (AHP) 
method and trend analysis was performed using Mann-Kendall test. Results indicated that mean 
annual integrated UV-VIS 200-440 radiation of north and northwest of the country has 
experienced a positive trend during the last four decades, but no significant trend was detected for 
other parts of the country. The study of the TCO also showed that in all eight clusters, the TCO 
trend during 1979 to 2000 has experienced a significant decreasing trend, but has increased from 
2001 onward. The linear regression equations between UV-VIS200-440 with TCO and CF were 
also constructed for all time scales (monthly, seasonal, annual) and the corresponding coefficients 
of correlation were calculated accordingly. The maximum annual correlation (r) was found 
between UV-VIS200-440 and the cloudiness (0.81). However, the highest r between UV-VIS200- 
440 and the TCO was 0.25. This result highlights that the variation in cloudiness plays more 
important role on determination of surface UV radiation than the changes in ozone concentration. 
This study did not investigated the influence of aerosol optical depth on UV radiation. This issue 
demands further investigation for evaluating the impact of aerosols on incoming surface UV 
radiation reaching the ground. 
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جغرافیا و مخاطرات محیطی. سال سیزدهم. شمارۂ چهل و نهم بهار ۱۶۰۳ صص ۲۶۳-۲۹۶ 
مقاله پژوهشی 
واکاوی چنددهه‌ای تغیبرات GRU‏ تجمعی روزانه فرابنفش سطحی در ایران و همبستگی آن با تغییرات 
ازون کلی و ابرناکی 
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حکیده 


تابش فرابنفش خورشیدی رسیده به سطح زمین تأثیر بیولوژیکی قابل‌توجهی بر 
سلامت انسان» حیوانات» گیاهان و اکوسیستم‌های زمینی و آبی دارد. پایش مستمر 
تغییرات شدت تابش فرابنفش و ازون IS‏ در ایران فقط در دو ایستگاه انجام می‌شود 
که در دوره‌هایی به خاطر مشکلاتی ازجمله کالیبراسیون و هزینه سنگین نگهداری با 
مشکلاتی مواجه شده است. این مشکلات انگیزه استفاده از دیگر منابع اطلاعاتی نظیر 
ماهواره‌ها و داده‌های باز تحلیل را در طیف‌های مهم الکترومغناطیسی خورشید فراهم 
می‌نماید. در پژوهش حاضر تابش تجمعی روزانه فرابنفش-مرئی (انتگرال کلی تابش 
فرابنفش-مرئی از طول‌موجح ۲۰۰ تا 460 نانومتر) حاصل از پایگاه داده بازتحلیل 
ERAS‏ ازون کلی و درصد ابرناکی مربوط به دوره آماری ۲۰۲۰-۱۹۷۹ در گستره 
ایران موردبررسی قرار گرفت. منطقه موردمطالعه با استفاده از روش خوشه‌بندی 
سلسله‌مراتبی به ۸ خوشه تقسیم و بررسی روند در سری داده‌ها با استفاده از آزمون 
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(UV-VIS200-440)‏ نشان داد در مقیاس سالانه در نواحی شمال و شمال غرب کشور 
(خوشه‌های ۵ ۰ ۷ و (A‏ روند مثبت وجود eagla‏ اما در ple‏ مناطق روند معنی‌داری 
مشاهده نشد. مطالعه سری زمانی ازون کلی کشور نیز نشان داد در هر ۸ خوشه مطالعه 
شده ازون کلی در سال‌های ۱۹۷۹ الی ۲۰۰۰ روند کاهشی معنی‌دار و از سال ۲۰۰۱ 
میلادی به بعد روند افزایشی معنی‌دار را تجربه نموده است. در این پژوهش به‌منظور 
تخمین تابش در نقاط فاقد ایستگاه» روابط رگرسیونی خطی بین تابش تجمعی(۲۰۰ 
الی ٤٤١‏ نانومتر) با ازون کلی و ابرناکی در مقیاس‌های ماهانه. فصلی و سالانه 
استخراج شدند. بررسی ضرایب همبستگی بین سری زمانی تابش فرابنفش-مرئی (۲۰۰ 
الی 8۶۰ انومتر) با ابرثاکی و ازون IS‏ نشان داد در هر سه مقیاس زمانی ماهانه 
فصلی و سالانه. GRU‏ با ابرناکی همبستگی قوی‌تری نسبت به ازون کلی دارد. به- 
طوری‌که در مقیاس سالانه بیشترین ضریب همبستگی تابش فرابنفش-مرئی با ابرناکی 
۱ و برای ازون IS‏ ۰/۲۵ به دست آمد. 

کلیدواژه‌ها: خوشه‌بندی سلسله مراتبی» پهنه‌بندی تابش روزانه» آزمون من-کندال» 


پایگاه داده ERAS‏ روند تابش فرابنفش. 
۱- مقدمه 


تابش اشعه فرابتفش خورشیدی (UV)‏ که به سطح زمین می‌رسد تأثیر قابل‌توجهی بر سلامت انسان؛ 
بهره‌وری کشاورزی, اکوسیستم‌های زمینی و آبی و کیفیت هوا دارد (داگلاس و همکاران ۲۰۱۱ "و ویلیامسون 
کارا اور eee leat‏ یک E‏ ات که تن ان ela cat‏ شیر Gedy pV:‏ وا 
فیلتر می‌کند. ازون استراتوسفری به‌طور موثر بخش زیادی از تابش فرابنفش خورشیدی را در طول‌موج‌های 
۰ تا ۳۱۵ نانومتر جذب می‌کند. که در غیر این صورت به‌طور بالقوه به حیات بیوسفر آسیب می‌رساند. در 
El‏ ند ۱۹۷۰ a‏ که اتسار مواد eal a‏ کل ما کسی بش ی رالد اعت کرو 
a‏ ازن شود. پروتکل مونترال در سال ۱۹۸۹ به اجرا درآمد و مواد متعددی را که باعث تخریب لایه ازون 
می‌شوند ممنوع کرد. از سال ۱۹۹۷ تا اواسط دهه ۲۰۰۰ آشکار شد که کاهش در ستون ازون (TCO) JS‏ 
تقریباً در تمام عرض‌های جغرافیایی غیرقطبی متوقف شده است (سازمان جهانی هواشناسی" ۲۰۰۷). 
ell‏ حال. سرعت بهبود تحت تأثیر تغییرات cleo‏ گردش جو و چرخه‌های تشکیل و تخریب ازون توسط 


نیتروژن و هیدروژن است. Ol jes‏ تابش فرابنفش خورشیدی دریافت شده در هر مکان خاص روی سطح زمین 
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علاوه بر میزان ازون در جو به موقعیت خورشید در افق و ابرناکی نیز بستگی دارد. در ایران مطالعات زیادی 
در زمینه GLU‏ فرابنفش انجام نشده است. در ادامه به چند نمونه از این مطالعات اشاره می‌شود. موقری و 
خسروی در سال ۱۳۹۳ در پژوهشی به بررسی توزیع مکانی شاخص پرتو فرابنفش در گستره ایران پرداختند. 
در این پژوهش به‌منظور محاسبه میزان پرتو فرابنفش. ٩۶‏ ایستگاه موردبررسی قرار گرفته است. نتایج نشان 
داد در مقیاس سالانه» ٤‏ درصد از مساحت کشور از شاعص پرتو فرابنفش متوسط برخوردار است که برحی 
شهرهای حاشیه دریای خزر را در برمی‌گیرد. ۵۶ درصد از سطح کشور که نیمه شمالی کشور را شامل می- 
شود. شاخص زیاد پرتو فرابنفش دریافت می‌کند و 4۲ درصد از مساحت کشور که عموماً استان‌های جنوبی 
کشور را در برمی‌گیرد. از میزان شاخص خیلی زیاد برخوردار است (موقری و خسروی. ۱۳۹۳). سبزی‌پرور 
و سیف‌زاده در سال ۱۳۹۹ با استفاده از مدل انتقال تابش لایه‌ای TUVS‏ به برآورد تابش تجمعی روزانه 
UVA‏ در سه حالت شرایط آسمان کاملاً صاف SUIS‏ ابری و واقعی در مناطق مرکزی ایران پرداختند. نتایج 
این پژوهش نشان داد بیشینه تابش UVA‏ در نیمه جنوبی منطقه موردمطالعه متمرکز است که در فصول گرم 
سال WE‏ در شرق و در فصول سرد سال در نواحی مرکزی و جنوب غرب منطقه مستقر است. در تحقیق 
دیگری سبزی‌پرور و سیف‌زاده (۱۳۹۹) به واکاوی تأثیر ابرناکی بر تابش فرابنفش تجمعی روزانه UVB‏ در 
مناطق Kis‏ و نیمه‌حشک dlp)‏ پرداختند. برآوردها در شرایط آسمان کاملاً صاف و ابری حکایت از حضور 
پیشینه تابش UVB‏ در فصل تابستان در نواحی جنوبی و جنوب شرقی منطقه موردمطالعه داشت. بر اساس 
نقشه‌های به دست آمده در این پژوهش. کمینه polis‏ روزانه UVB‏ نیز در تمامی شرایط موردبررسی در 
فصل زمستان مشاهده شد. سبزی‌پرور و همکاران (1۹۹۹) با استفاده از مدل تابش مالتی استریم (چند چشمه- 
ای) و جو چندلایه‌ای به ارائه میانگین‌های درازمدت تابش‌های بیولوژیکی UVB UVA‏ آفتاب‌سوختگی و 
آب‌مروارید چشمی در مقیاس ۱/۲۵۴« ۱/۲۵۴ اقدام نمودند. در این تحقیق نشان داده شده که در توزیع 
جغرافیایی تابش‌های بیولوژیکی موقعیت سمت‌الرس خورشید (زاویه پرتوهای خورشید با خط عمود بر 
محل) نقشی مهم‌تر از ازون ایفا می‌کند. فرج‌زاده اصل و همکاران (۱۳۹۳) به تحلیل تغییرات تابش فرابنفش 
در چهار مقیاس زمانی روزانه. ماهانه» فصلی و سللانه در منطقه اصفهان پرداختند. نتایج حاصل از برازش 
روند تغییرات فصلی نشان داد مقدار تابش در فصل بهار روندی افزایشی داشت. درحالی‌که این روند برای 
دیگر فصل‌ها کاهشی بود. همچنین نتایج نشان داد بیشترین شاخص تابش فرابنفش (UV Index)‏ به مقدار 
۵ در اواسط تابستان و کمترین Ol‏ به مقدار ۰/۵ در زمستان رخ داد. رستم‌پور و همکاران (۱۳۹۱) در 
پژوهشی به بررسی میزان شدت تابش ماهانه UV-A‏ در شهر همدان پرداختند. این مطالعه به‌صورت مقطعی 
و در مدت یک سال انجام شد. نتایج این پژوهش نشان داد بيشینه میانگین UV-A‏ در شهریور با میزان 
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کمینه میانگین تابش UV-A‏ نیز حدود ۱۸/ از این مقدار مجاز تجاوز کرد. علت این اختلاف به ارتفاع زیاد 
شهر همدان از سطح دریا نسبت داده شده است. 

پایش تغییرات شدت GLU‏ فرابنفش و ازون در Ol gl‏ به‌صورت محدود فقط در دو ایستگاه انجام می‌شود 
که به‌صورت غير پیوسته و نامنظم است. مطمئناً در کشور پهناوری همچون ایران» اندازه‌گیری نامنظم تنها در 
دو ایستگاه نمی‌تواند کافی باشد و این امر ضرورت استفاده از داده‌های شبکه‌ای جهت بررسی میزان تابش 
فرابنفش و ازون» همچنین بررسی تغییرات این تابش به‌عصوص در سال‌های اخیر را نشان می‌دهد. داده‌های 
5 نسل پنجم داده‌های بازتحلیل مرکز پیش‌بینی‌های میان‌مدت اروپا ECMWF‏ از اقلیم جهان برای ۶ تا 
۷ دهه اخیر است. ERAS‏ تخمین‌های ساعتی را برای تعداد زیادی از مقادیر جوی, امواج اقیانوسی و سطح 
زمین فراهم می‌کند. برآورد عدم اطمینان توسط روش Ensemble‏ تحت ده بخش در فواصل سه‌ساعته انجام 
می‌شود. این قبیل برآوردهای عدم اطمینان از نزدیک با محتوای اطلاعاتی سیستم مشاهده موجود مرتبط است 
که با گذشت زمان به‌طور چشمگیری تکامل یافته است. چند نمونه از پژوهش‌هایی که با استفاده از داده‌های 
تابش فرابنفش ERAS‏ انجام شده است در ادامه ذکر می‌شود. شیا و همکاران (۲۰۲۱) به بررسی اثر افزایش 
ابر ناشی از تخریب ازون بر میزان UV‏ رسیده به سطح در سیبری در فصل بهار با استفاده از داده‌های ERAS‏ 
پرداختند. آن‌ها نتیجه گرفتند که از بین رفتن لایه ازن استراتوسفر باعث افزایش اشعه فرابنفش رسیده به سطح 
می‌شود. در مقابل افزایش ابرهای مرتفع ناشی از تخریب ازون منجر به کاهش UV GRE‏ سطح می‌شود. آن- 
ها با یک تجزیه‌وتحلیل ترکیبی به این نتیجه رسیدند که اثر افزايش ابر قوی‌تر از اثر از بین رفتن ازون 
استراتوسفر در منطقه موردمطالعه است. راکسسات و همکاران (۲۰۲۱) به پیش‌بینی دقیق تابش فرابنفش 
سطحی برای کاربردهای بالینی با شبکه عصبی در استان ناکون پاتوم تایلند پرداختند و دقت مدل‌ها را با 
استفاده از داده‌های تابش فرابنفش, ابرناکی و ازون کلی ERA5‏ مورد ارزیابی قرار دادند. 

در زمینه بررسی روند طولانی‌مدت داده‌های فرابنفش نیز مطالعاتی انجام شده است. گلندورف و همکاران 
(۲۰۰۶)" به بررسی روند تابش فرابنفش در سری زمانی دو ایستگاه در یونان و فنلاند پرداختند. نتایج نشان 
داد روند معنی‌داری در این دو سری زمانی در داده‌های تابش فرابنفش مشاهده نمی‌شود. سمدلی و همکاران 
(۲۰۱۲)* به بررسی روند تغییرات تابش فرابنفش سطحی در انگلستان پرداخت. نتایج این پژوهش نشان داد 
از سال ۱۹۹۳ تا ۲۰۰۸ روند قابل‌توجهی در میزان حداکثر شاخص UV‏ روزانه AVI)‏ در هر دهه) دیده می- 


شود. اما هیچ ارتباطی بین تغییرات UV‏ سطح و مقادیر ازون کل وجود نداشت که نشان می‌دهد میزان UV‏ 
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بیشتر تحت تأثیر تغییرات در پوشش ابر است. برنارد و استایرل (۲۰۲۰) به بررسی روند در سری زمانی 
تابش فرابنفش ۲۰۱۸-۱۹۹۰ در قطب جنوب پرداختند. نتایج این مطالعه شواهدی ارائه کرد که نشان می‌دهد 
شاخص (UVD‏ اندازه‌گیری شده توسط طیف‌سنج‌های 517۷-100 در سه سایت قطب جنوب (قطب جنوب؛ 
ارتفاعات Arrival‏ و ایستگاه پالمر) در حال کاهش است. به‌عنوان‌مثال» یک روند نزولی قابل‌توجه ۵/۵ 
درصد در هر دهه (۹۵/ سطح اطمینان) در ارتفاعات Arrival‏ برای زمستان (دسامبر تا فوریه) گزارش شد. 
فانتولکیس و همکاران (۲)۲۰۱۸ با استفاده از داده‌های بیست‌وپنج سال اندازه‌گیری طیفی UV-B‏ در کاناد؛ 
ارویا و ژاپن» به بررسی روند و اثرات تغییرات در ازون» آئروسول, ابرها و آلبیدو بر میزان تابش UV-B‏ 
پرداختند. نتایج پژوهش آن‌ها OLS‏ داد در نیمکره شمالی تغییرات طولانی‌مدت در تابش UV-B‏ رسیده به 
سطح. در مکان‌های مختلف به‌طور قابل‌توجهی متفاوت است. محرک‌های hel‏ این تغییرات؛ تغییرات در 
آئروسول و ازون IS‏ گزارش شد. فانتولکیس و همکاران (۲۰۲۰)" در پژوهشی به بررسی روند در سری 
زمانی تابش فرابنفش در اروپا پرداختند. نتایج تجزیه‌وتحلیل ازون کلی و داده‌های طیفی 0۷ ثبت شده در 
چهار ایستگاه ارویایی طی شال‌های: ۱۹۹۰ تا ۲۶۱۷ نشان داد تغییرات درازمدت در UV‏ عمدتاً ناشی از 
تغییرات در آئروسول» ابرناکی و آلبیدو سطح است. درحالی که تغییرات در ازون کلی نقش کمتری دارد. اون 
و همکاران (۲۰۱۹)* به بررسی ویژگی‌های روزانه» فصلی و سالانه UVB GSR‏ (۳۱۵-۲۸۰ نانومتر) و 
UVA‏ (1۰۰-۳۱۵ نانومتر) ازون IS‏ و غلظت آئروسل‌ها از سال‌های ۲۰۱-۲۰۰6 پرداختند. نتایج پژوهش 
Lol‏ نشان داد هیچ روند قابل‌توجهی در GLU‏ فرابتفش برای دوره موردمطالعه در مقیاس سالانه و فصلی 
وجود نداشت. همین امر در مورد عوامل موثر بر تابش فرابنفش, ازون کلی. GSR‏ و عمق نوری آثروسل 
صادق بود. 
تحقیقاتی که تا کنون در کشور انجام شده است. به علت نبود ایستگاه‌های اندازه‌گیری زمینی تغییرات 
بلندمدت میزان تابش فرابنفش و وجود روند احتمالی در آن موردبررسی قرار نگرفته است. در تحقیق حاضر 
با استفاده از داده‌های بازتحلیل ERAS‏ اقدام به پهنه‌بندی تابش فرابنفش (انتگرال کلی GLU‏ فرابنفش از 
طول‌موج ۲۰۰ تا ٤٤0‏ نانومتر که در کل مقاله به خاطر تلخیص تابش فرابنفش نوشته می‌شود) و ازون کلی و 
ابرناکی و همچنین بررسی تغییرات میزان تابش فرابنفش و ازون IS‏ و ابرناکی طی سال‌های ۱۹۷۹ تا ۲۰۲۰ 
در گستره Ul pl‏ می‌شود. ذکر این نکته ضروری است که در این تحقیق امکان جداسازی تابش ناشی از ۶۰۰ تا 
۰ نانومتر در داده‌های بازتحلیل ERAS‏ وجود نداشت و به این دلیل محدوده طبف موردمطالعه (11۰-۲۰۰ 
نانومتر) فرابنفش نام گذاری شد. 
Bernhard and Stierle, 2020‏ 1 
Fountoulakis et al, 2018‏ 2 


3 Fountoulakis et al, 2020 
4 Aunet al, 2019 
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-Y‏ منطقه موردمطالعه 


کشور ایران با وسعت ۱۹۶۸۰۰۰ کیلومترمربع در نیم‌کره شمالی» در قاره آسیا و در قسمت غربی فلات 
ایران واقع شده است. این کشور. بین دو نصف‌النهار "44 و 14۳ شرقی و دو مدار ۲۵۳ و ۶۰۳ شمالی قرار 
گرفته eal‏ حدود ٩۰‏ / خاک ایران در محدوده فلات ایران واقع شده انیت و کشوری کوه ‌ستانی 


محسوب می‌شود. علاوه بر گسترش عرضی» رشته‌کوه‌های البرز و زاگرس مهم‌ترین BE‏ را در آب‌وهوای 
ایران دارند. افزون بر این» دو گستره آبی بزرگ در شمال و جنوب. رژیم‌های دمایی و بارشی خاصی را در 
کشور ایجاد می‌کند. به‌طوری که شرایط رطوبتی در این سرزمین از زمین‌های فوق‌العاده خشک تا زمین‌های 
بسیار مرطوب را شامل می‌شود. در این پژوهش» کشور ایران و نواحی اطراف OF‏ به‌صورت شبکه‌ای با فاصله 
نقاط ۲۵ کیلومتر موردمطالعه قرار گرفت. تعداد نقاط موردبررسی در این پژوهش 1۵4۳ نقطه است که توزیع 
نقاط شبکه در شکل ۱ نشان داده شده است. 


50°0'0"E 55°0'0"E 60°0'0"E 


I A هه‎ sid t 1 7 
FOO IOI IT ۱ FOI III III III III III OI OI IO IOS ISAAK 
FOI KKK ORK KOK OK OKO KHKKKAKHKKKKKKKKHKHNHHHHKHKKKKKKHNKN KKK 
IOI OOO II I IK FOI III OI III IOI I IOI IOI IO III 
SRO OOO I tor k ۳ E E E با‎ IFA 
FOI III III E EEEE EEEE 
0 چا هچ ری ی هد‎ IIIT IOI 
مق ال دز‎ k STII وی‎ 
FO I I IIIT TOR TO KR IO ** 5 3 
FORO OTTO OTTO RI IE 
۳ : 
wee EE EEE EEE EE EEE EEE EE EEE ee ee ee 
FO OI III III IOI III III IIIT III Ir 
FOI چا چا مج هه مر چم‎ III چا با‎ IOS IIIA ISIS AIA ISIS IAI AISI SSSI I IIIS SI ISA IIIA ASI IIASA 
FOO را‎ 
FOR III IIIS III III IOI IDI IOI II III IOI IOI IOI IOI IOI IO IOI OI IIIA III AIA OS IAA ۷ ۷ 
Aaaa با ار‎ IDIOTS IIIT II IOI III 
FOI III IOI چا چا با‎ IIIS III III III III TOI III 
هچ بر 6 ۳ ار‎ III II III III IOI III III III IO IOI III II tt 
ها جر هچ و چا له‎ OK OK OK OKC چا‎ Û DK KOK DK OK DK OK OK OK OK DK DK KK DK OK Û Û Û OK DK OK O II IOI IOI TOO TOD IORI A SOTA DSI SE IOI OB AAAI 
FORO IISA OSS IIA III III I IAA ASIII III SII IAI ALI I AISA AS AIA IA A. 
FOI III III III III IOI DK KOK OOK I IORI OK O DÛ TOK IIA IOS IA II IAAI IIIS AISI E SIA ISI IA IAA. 
FIO III IIIS OI OTS I IOS IOI I IIA IIA III ISIS ASI ISI IAI III AISA SI AIA I OIA IA IAA. 
FI III III II IOI III IOI IOI OOS I IOS IOI II IIA IIA IAS III SAI III III SII III AA IAAI SI IAI AA. 
FO یی و و‎ III IOI II 
FOI چم‎ OK OK OK OK KKK OK OK OK OK OOK KOK OK DK Û OK چم ها‎ III OK Û O KO KO OK OK OK IOI OK OK KOK OK OK OK OI TOI II: 
FOI III III IOI ICICI IOI IOI IOI OI با‎ IA IIA III I IA IIIA IIS IAS ASAI II ISI با‎ ASAI I IASC AI IAS ACS ۷ ۰ ۷ ۸ 
FIR I IOI IOI I IO IOI IO IOS IA II 
FIT IOS IIIA IIA AA IIA SIA ABI IA BIIIA SAS AII IIIS IIA IIA I ASIA ANIA ISA A IA 
FOI III یی ویب ما‎ OI OI III 
FORO III III OI III OI IOI IOI OOS I IOS AOI OSI IIA I IIA IAI IIIA ASIII III III IAI AAA IIA AAI IAI OAA. 
اج اج ها چا ی سر‎ IIGIIDISIDISIIOIIOIIDIDIDIDIUIOIIOIDIDIDIDIDIDIOIIOIIDIDIDIDIIOIIDIDIDIDIDIOIDIO IOI IOI IGIDISI I TOI KKK 
FOI هب هچ ی‎ III چا چا‎ OK OK OK OK OK KOK OK OK چا‎ IOI III IOI IOI IOI OI IOI OIA OT AIA AIO AIS SIA IAA. 
FI III III III IO IOI IOI IO III TOSI EEEE OS ISAS IA IIIA IIIA AI ISI IAI IIA II ASAI ISAACS IAAI AA 
FIO TOI II III III OI OK OK I IOS IOI SI SIA IIIA IIIA II IIA I ISI SII III IIIA IISA ISAS AISA. 
FO IORI IOS ITOK IOI OKI ORR IOI IOI IO IIIT IA IIA II II SII IA III I IIIA I IAI IAI SI جر‎ IAA. 
FOI IO IOI IR FIO III III III III II ORI IIIA IOI IIS ISI ISIS SSIS AISA IS AN 
FOI IOI I RIO III III III III III IO III OI OIA TIAA IH 
FOI Ok 1 ماس ی و‎ III KKK 
۱ TIO Ik ۳ III IOI OI IOI IIA III IAAI OS III IAAI IA. 
FOI IOI KK k ETETETT III III III III IOI III I IOI OIA IOI A IORI III IIIA 
FOI ITI TI E E OI IOI IIIS OK IORI IOI OK IOI II III ISI SIA E E IIA ISIN. 
1۱ ttt E E E E E E E E E E E E E E E E E E E E E E E E E E IOI Itt 
ETETETT EE EEEE EEE & 
FOI Rk FOI III IOI III IOI IO IOI IO IOI OITA II AIA SA SIDI 
FOI III III FOI III III III III IOI OA ITA SOA IIA INA 
FO II III IOI III III k EEEE EEE E A N A EEE EE A E AA EEEE EEE A A A A A E 
مج با مب سا‎ III IO - و‎ IOI IOI E EEEE IOI IO IOI OI IIA IIA 
FO III IIIT II OTTO FEN 
CIT Tt i i tk ttt tlie 
FOO III OI IO RIO 
EE Ito 
و‎ IIIT IOI II 
FOI III OI IOI IOI IK RR RH 
۱ k k 7 KKH 
FORO Rk ane i ea aad 
FOI III k ttrt 


50°0'0"E 55°0'0"E 60°0'0"E 


LJ Kilometers 
0 75150 300 


۳- مواد و روش‌ها 

ERAS داده‌های‎ 

پایگاه اطلاعاتی ECMWF‏ یکی از مهم‌ترین پایگاه‌های داده بازتحلیل در جهان است که موردتوجه 
بسیاری از محققین قرار دارد (دی و همکاران» jl OS‏ تحلیل داده‌های مدل را با استفاده از قوانین فیزیک 
با داده‌های مشاهده‌ای از سراسر جهان ترکیب می کند و یک مجموعه داده کامل و سازگار جهانی ایجاد می- 
کند. این اصل که داده‌گواری " نامیده می‌شود مبتنی بر روشی است که توسط مراکز عددی پیش‌بینی آب‌وهوا 
استفاده می‌شود. که در آن هر چند ساعت (۱۲ ساعت در ECMWF‏ پیش‌بینی قبلی با مشاهدات جدید در 
دسترس» با روشی بهینه ترکیب می‌شود تا بهترین تخمین جدید از وضعیت جو تولید شود. ERAS‏ نسل 
پنجم بازتحلیل ECMWF‏ برای آب‌وهوا و اقلیم Slee‏ در ٤‏ تا ۷ دهه گذشته است. در این پژوهش داده‌های 
تابش فرابنفش بازتحلیل ERAS‏ ابتدا به‌صورت میانگین ماهانه و با رزولوشن مکانی ۲۵ VOR‏ کیلومتر برای 
گستره ایران از آدرس https://cds.climate.copernicus.eu‏ دریافت شد. این داده‌ها مربوط به انتگرال کلی 
تابش فرابنفش و بخشی از تابش visible‏ از طول‌موج ۲۰۰ تا 460 نانومتر می‌باشند که از سال ۱۹۷۹ تا ۲۰۲۰ 
JS)‏ دوره موجود) دانلود گردید. با توجه به اینکه محدوده طول‌موج‌های visible‏ از طول‌موج‌های مربوط به 
تابش فرابنفش قابل‌تفکیک نیست. تابش فرابنفش (UV-VIS200.440)‏ در مقاله به کل این محدوده از ۲۰۰ تا 
۰ نانومتر اشاره دارد. همچنین داده‌های ازون و ابرناکی به‌صورت میانگین ماهانه و با رزولوشن مکانی ۲۵ 
# کیلومتر دانلود و پس از کنترل کیفی با استفاده از نرم‌افزار R‏ و اطمینان از نرمال بودن و عدم وجود داده- 
Gh‏ پرت مورداستفاده قرار گرفت. این داده‌ها توسط نرم‌افزار Open Netedf‏ به‌صورت شبکه‌ای برای 101۳ 
نقطه برای گستره کشور استخراج گردید. 

خوشه‌بندی سلسله مراتبی 

alls ed‏ ی ای yy‏ است که به کمک وه گ ها با صجفت‌ها ی مشاهانی: آن wily‏ گروه‌های 
مشابه طبقه‌بندی می کند. انتخاب ویژگی‌های مناسب برای این کار یکی از مسائل مهمی است که باید در نظر 
گرفته شود. از طرف دیگر استانداردسازی داده‌ها نیز مطرح است تا مقیاس اندازه‌گیری صفت یا ویژگی‌ها 
باعث انحراف تابع فاصله نشود. در این پژوهش نرمال بودن داده‌ها با استفاده از نرم‌افزار R‏ آزمون شد و برای 
انجام روش خوشه‌بندی سلسله مراتبی از نرم‌افزار SPSS‏ استفاده شد و خوشه‌بندی با استفاده از داده‌های 


میانگین تابش فرابنفش سالانه» میانگین ابرناکی سالانه و میانگین ازون سالانه و عرض جغرافیایی و ارتفاع از 


1 Dee et 1 

2 Data assimilation 
3 Features 

4 Attributes 
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سطح دریا انجام شد. درروش خوشه‌بندی سلسله مراتبی از فاصله اقلیدسی جهت محاسبه فاصله استفاده 
گردید (رابطه ۱). هر چه همبستگی بیشتر باشد فاصله کمتر در نظر گرفته می‌شود. به‌طوری‌که با اندازه‌گیری 
فاصله اقلیدسی متغیرهای میانگین تابش فرابنفش سالانه» میانگین ابرناکی سالانه و میانگین ازون سالانه و 
عرض جغرافیایی و ارتفاع از سطح دریا مربوط به نقاط موردمطالعه» همبستگی و همانندی نقاط مشخص 
می شود. 
Deuce = C= 0)‏ 
شکل ۲ مربوط به ag‏ درون‌یابی شده L)‏ روش IDW‏ مربوط به میانگین سالانه ابرناکی (راست) و 
میانگین سالانه ازون کلی (چپ) است. شکل ۳ پهنه تابش تجمعی روزانه 200-440 UV-VIS‏ را نشان می‌دهد. 
درنهایت منطقه موردمطالعه به ۸ خوشه تقسیم شد که پهنه‌بندی خوشه‌های تعیین شده در شکل ٤‏ آورده شده 


است. در این پژوهش پهنه‌ها با استفاده از نرم‌افزار ARC GIS‏ ترسیم شد. 
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شکل ۲- توزیع جغرافیایی میانگین سالانه ازون کلی و ابرناکی (۱۹۷۹-۲۰۲۰) منطقه موردمطالعه با استفاده از 


روش 10۷۷ 


سال سیزدهم واکاوی چنددهه‌ای تغییرات تابش تجمعی روزانه فرابنفش سطحی در ایران و... Yo\‏ 
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شکل ۳- توزیع حغرافیایی میانگین سالانه تابش فرابنفش تحمعی روزانه (۱۹۷۹-۲۰۲۰) منطقه موردمطالعه با 
استفاده از روش IDW‏ 
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شکل ag -٤‏ خوشه‌بندی سلسله مراتبی منطقه موردمطالعه با استفاده از روش IDW‏ 
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آزمون من -کندال 
جهت بررسی وجود روند در سری زمانی داده‌های GLU‏ فرابنفش» ازون و ابرناکی از آزمون من-کندال در 
بستر نرم‌افزار آماری R‏ استفاده شد. آزمون من -کندال» آزمون ناپارامتریک رتبه مبنا جهت بررسی روند است. 
در این آزمون تحت فرض Ho‏ داده‌ها از سری که مستقل و دارای توزیع یکسان می‌باشند گرفته شده‌اند. 
آزمون من -کندال با استفاده از روابط ۲ و ۳ انجام می‌شود : 
)۲( 0 < ,1+ 
sgn(x; — x;), where sgn(x) =} 0,x = 0‏ مین s = LT‏ 
—1,x< 0‏ 
که واریانس فرآیند از معادله (۳) تعیین می‌شود. 
o2 = —[N(N - 1)(2N +5) - DM tit; - 1)(2t + 5)] e‏ 


که در cpl‏ رابطه m‏ تعداد دسته‌ها با داده‌های گره‌دار و ti‏ تعداد داده‌های گره‌دار در هر دسته M‏ می‌باشد. 

Sen’s slope روش‎ 

OAT) Sen‏ با توسعه و بسط یک سری مطالعات آماری که Thiel‏ )+140( به انجام رسانده بود یک 
روش ناپارامتری را جهت تحلیل سری‌های زمانی ارائه نمود. این روش نیز همانند روش من- کندال از 
تحلیل تفاوت بین مشاهدات یک سری زمانی بهره می گیرد. اساس این روش بر محاسبه یک شیب میانه می- 
باشد. مراحل LIS‏ انجام این آزمون به شرح زير می‌باشد: محاسبه شیب بین هر جفت osla‏ مشاهده‌ای با 
استفاده از رابطه :٤‏ 


Q =% (é) 


ts 
s و , × به ترتیب داده‌های مشاهده‌ای در زمان‌های ] و هو ] یک واحد زمانی بعد از زمان‎ 2, wi که در‎ 

EE اعمال واک ا ای کر دنت کاده اها ایک سر رما از کدی‎ Lal 

حاصل می‌آید که از محاسبه میانه این سری زمانی شیب خط روند (Q med)‏ به دست می‌آید. مقدار مثبت Q‏ 

Smead‏ از صعودی بودن روند و مقدار منفی آن نشان دهنده نزولی بودن روند است. 

در این پژوهش محاسبه ضرایب همبستگی و استخراج معادلات رگرسیون چند متغیره با استفاده از نرم‌افزار 

SPSS‏ انجام شد. 


شال سیزدهم واکاوی جنددهه‌ای تغییرات تابش تجمعی روزانه فرابنفش سطحی در ایران و... yor‏ 


-٤‏ نتایج و بحث 


۱-۶ مقیاس فصلی 


وا کے عفر ات شالانه Salo Ge‏ رو ۱۵3۷ ارین کل ی SUA‏ در طرل 99 bj‏ 
آماری LY‏ ساله» در ۸ خوشه در محدوده موردمطالعه بررسی شد و وجود یا عدم وجود روند با استفاده از 
آزمون من -کندال مورد آزمون قرار گرفت. با توجه به شکل ۵ در فصل زمستان تغییرات میزان ابرناکی در 
تمام خوشه‌های موردبررسی منفی بوده است اما این روند کاهشی به‌جز در خوشه‌های ۱ و ۸ در سایر 
خوشه‌ها معنی‌دار بوده است. در رابطه با ازون کلی و 0۷-۷180440 روند معنی‌داری در هيچ‌کدام از 
خوشه‌ها در فصل زمستان وجود نداشت. تغییرات میزان ابرناکی در فصل بهار نیز در تمام خوشه‌ها کاهشی 
ات اما فقظ دو شر نها Gg Oe oY‏ ای رون Sale‏ مان رده اس رات و55 ۱0۷2۲ 
در فصل بهار در تمام خوشه‌ها افزایشی است اما فقط در خوشه ۵ این روند معنی‌دار است و در فصل بهار 
هیچ‌گونه روند معنی‌داری در رابطه با میزان ازون کلی وجود ندارد. تنها روند معنی‌دار موجود در میزان 
ابرناکی در فصل تابستان مربوط به خوشه ۱ است که می‌تواند مربوط به مناطق جنوب شرقی ایران و پدیده 
igh‏ رات ای تاک هو وهای Vy lls‏ وی MAIS‏ اش انا ورتم دای سیگ 
A es EVES ies‏ ون هر فیهای تم ۷و دازا ری آفرانشی مک وان pica‏ 
میزان ازون IS‏ در فصل تابستان در تمام خوشه‌ها روند افزایشی معنی‌دار داشته است. همچنین روند ازون 
کلی در فصل پاییز در تمام خوشه‌ها به‌طور معنی‌دار افزایش یافته است. 0۷-۷1820044 در خوشه‌های ۸ ود 
روند افزایشی و در خوشه‌های ۱ و ۲ کاهشی org‏ است و در ply‏ خوشه‌ها روندی مشاهده نشد. بررسی 
داده‌های ابرناکی نشان می‌دهد که در فصل jul‏ میزان ابرناکی در خوشه‌های V‏ و ۸ کاهش adh‏ است. هر 
tse‏ این اهن تون قعش ذاری فار اما به ف زم غات EAN‏ ودن رونت gS VV IS gpg‏ این 
cll‏ ا اک تال هجیتن بقل رو اه یش ی شون ماه a‏ 
۲ در فصل پاییز روند افزایشی ازون IS‏ می‌باشد زیرا Gb‏ جدول ۳ در این دو خوشه روند افزایشی ازون 
کلی شیب بیشتری نسبت به ple‏ خوشه‌ها دارد. 

els‏ ضرایب همبستگی () تابش با ازون و ابرناکی (جدول )٤‏ نیز نشان داد ۲ مربوط به 0۷-۷1520440 و 
ابرناکی بسیار بیشتر از 0۷-۷152044 و ازون JS‏ است. درمجموع. بیشترین ضریب همبستگی مربوط به 
ابرناکی و ed ON VS‏ بیشترین ضریب همبستگی ابرناکی و وو ۷18 0۷ مربوط به فصل پاییز و 
ste‏ که اراق days‏ بالات راز اراک اتف اشد 5 کر کو یب هتکن UV-VIS sian 5 Stal)‏ 


yog‏ جغرافیا و مخاطرات محیطی شمارة اول 


خورشید در فصل زمستان باشد. بیشترین همبستگی ازون با تابش فرابنفش مربوط به فصل تابستان در خوشه 
1 است که به دلیل ابرناکی بسیار کم در این فصل در مناطق کوبری میزان تابش فرابنفش نسبت به فصول 
دیگر بیشتر تحت تأثیر تغییرات ازون و زوایای سمت‌الرآس کمتر بوده است. کمترین ضرایب همبستگی ازون 
و 0۷-۷180 مربوط به فصل زمستان و پاییز است. این نتایج با نتایج دیو پریز "و همکاران نیز که در 
سال ۲۰۱۹ به بررسی ارتباط بین تابش فرابنفش خورشیدی (UV)‏ و ازون پرداختند هم سو است. آن‌ها نتیجه 


گرفتند که بیشترین ضریب 8 بین ازون و تابش فرابنفش مربوط به شرایط آسمان صاف به مقدار ۰/٤۵‏ بود. 


ا soore SS'Q0E‏ شتا نا 


E 
z 
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شکل 0- روند فصلی تابش فرابنفش روزانه )390-440 (UV-VIS‏ ازون کلی و ابرناکی در خوشه‌های موردمطالعه 


1 Du Preez 


YOO 


واکاوی جنددهه‌ای تغییرات تابش تجمعی روزانه فرابنفش سطحی در ایران و... 


سال سیزدهم 


جدول ۱- آماره سطح معنی‌داری (آزمون من -کندال) و شیب Sem‏ میانگین‌های فصلی تابش فرابنفش (۱۹۷۹-۲۰۲۰) 


Winter Spring Summer Autumn 
Sen's slope p- Sen's slope p- Sen's slope Sen's slope 
Cluster\Test p-value Kj/m?.year value Kj/m’.year | value Kj/m’.year | p-value Kj/m’.year 
< 

Cluster 1 0.530 0.444 0.696 0.200 0.747 -0.045 0.0001 -0.906 
Cluster 2 0.135 0.696 0.053 0.852 0.016 0.182 0.001 -0.865 
Cluster 3 0.948 -0.073 0.222 0.975 0.474 0.164 0.095 -0.612 
Cluster 4 0.329 -0.646 0.233 1.127 0.075 0.707 0.588 -0.202 
Cluster 5 0.795 -0.138 0.022 1.402 0.046 1.007 0.802 0.093 
Cluster 6 1.000 -0.002 0.076 1.994 0.008 1.871 0.020 0.546 
Cluster 7 0.172 0.517 0.054 1.893 0.023 2.226 0.072 0.803 
Cluster 8 0.729 -0.315 0.146 1.931 0.044 2.860 0.009 1.282 


جدول -Y‏ آماره سطح معنی‌داری (آزمون من -کندال) و شیب Sem‏ (درصد در سال) میانگین‌های فصلی ابرناکی 


)۱۹۷۹-۲۰۲۰( 


Winter Spring Summer Autumn 

É 2 5 2 5 2 5 2 5 

i å 1 1 1 1 
Cluster 1 0.435 -0.1 0.258 -0.1 0.041 0.0 0.270 0.1 
Cluster 2 0.046 -0.1 0.014 -0.1 1.000 0.0 0.307 0.0 
Cluster 3 0.018 -0.2 0.032 -0.1 0.153 0.0 0.518 0.1 
Cluster 4 0.003 -0.3 0.076 -0.1 0.791 0.0 0.897 0.0 
Cluster 5 0.048 -0.2 0.022 -0.1 0.745 0.0 0.845 0.0 
Cluster 6 0.027 -0.2 0.114 -0.1 0.260 -0.1 0.269 -0.1 
Cluster 7 0.005 -0.2 0.042 -0.2 0.410 -0.1 0.179 -0.1 
Cluster 8 0.056 -0.1 0.091 -0.1 0.242 -0.1 0.072 -0.1 


۲ آماره نی‌داری (آزمون ۰ -کندال) و شیب روند 560 (دابسون در سال) میانگین‌های‎ eee 
oe معنی‎ 


ازون کلی (۱۹۷۹-۲۰۲۰) 


bs Winter Spring Summer Autumn 

Š a 20 é AT a E ۳ 20 

z E 3 = 3 و 3 = 3 = و‎ 

8 | 4 261 á اه‎ 4 2 E 8 
Cluster 1 0.752 0.097 0.280 0.169 < 0.0001 0.314 < 0.0001 0.244 
Cluster 2 0.986 0.001 0.583 0.094 < 0.0001 0.301 < 0.0001 0.253 
Cluster 3 0.965 -0.002 0.828 0.029 0.000 0.249 0.000 0.214 
Cluster 4 0.850 0.031 0.713 -0.034 0.002 0.216 0.002 0.199 
Cluster 5 0.816 0.052 0.423 -0.086 0.006 0.219 < 0.0001 0.216 
Cluster 6 0.853 0.019 0.696 -0.057 0.003 0.218 0.004 0.211 
Cluster 7 0.828 0.029 0.862 -0.025 0.003 0.211 0.002 0.201 
Cluster 8 0.931 0.013 0.948 -0.017 0.030 0.202 0.006 0.174 


ماه اون 
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جدول 4- ضرایب همبستگی فصلی بین تابش فرابنفش تجمعی روزانه با ازون کلی (TCO)‏ و ابرناکی (CF)‏ 


Clusters 


Cluster 1 
Cluster 2 
Cluster 3 
Cluster 4 
Cluster 5 
Cluster 6 
Cluster 7 
Cluster 8 


and CF UV-VIS200-440 


and TCO UV-VIS 200-440 


Winter 


Spring Summer Autumn 
-0.83 
-0.86 
-0.93 
-0.94 


-0.91 
-0.83 
-0.88 
-0.92 


Winter 

0.2 
0.03 
-0.13 
-0.12 
-0.23 
-0.31 
-0.34 
-0.28 


Spring Summer Autumn 
-0.37 0.24 -0.4 
-0.43 0.24 -0.43 
-0.48 0.13 -0.49 
-0.48 0.31 -0.44 
-0.44 0.41 -0.33 
-0.41 0.58 -0.19 
-0.37 0.52 -0.11 
-0.34 0.39 -0.05 


در جدول ۵ روابط رگرسیون خطی چند متغیره مربوط به خوشه‌های موردمطالعه در هر فصل ارائه شده 


کیلوژول بر مترمربع بر روز است. 
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واکاوی جنددهه‌ای تغییرات تابش تجمعی روزانه فرابنفش سطحی در ایران و... 


سال سیزدهم 


حدول 0— معادلات رگرسیون خطی جندمتغیره مربوط به تابش فرابنفش )300-440 (UV-VIS‏ ازون کلی TOC‏ و 


ابرناکی CF‏ در مقیاس فصلی 
Clusters Winter Spring‏ 

Cluster 1 1889.081+UV-VIS290.440=0.279TCO-679.4CF | yy Wisng ayg=-0.217TCO-854.13CF+3074.946 
UV-VIS 509-449=0.224TCO- 

Cluster 2 627.436CF+1738.858 | UV-VIS2009.449=-0.02TCO-909. 195CF+2953.504 
UV-VIS 590.449=0. 147TCO- 

Cluster 3 457.288CF+1602.044 | UV-VIS200-449=0.142TCO-1051.38CF+2910.157 

17۷-۷ 19200402۰ 1 82TCO- UV-VIS 200-449=0.4 1 2TCO- 

Cluster 4 399.716CF+1440.785 1285.659CF+2849.056 
UV-VIS 200-440=-0.006TCO- 

Cluster 5 407.444CF+1415.160 | UV-VIS290-449=0.117TCO-1178.85CF+2905.706 

UV-VIS 200-440=-0. 1 5 1 TCO- UV-VIS 200-449=0.2 1 2TCO- 

Cluster 6 416.739CF+1420.317 1213.679CF+2742.818 
UV-VIS 200-440=-0.560TCO- 

Cluster 7 545.662CF+1591.753 | UV-VIS200-440=0.64TCO-1280.567CF+2610.607 

UV-VIS 200-440=-0.8 1 6TCO- UV-VIS 200-449=0.723TCO- 

Cluster 8 711.382CF+1711.114 1628.137CF+2695.511 

Summer Autumn 

UV-VIS200-449=-0.779TCO- 

Cluster 1 UV-VIS300-449=0.27TCO-424.669CF+2944.653 487.195CF+2407.587 
U V-VIS200-440=0.47 1TCO- 

Cluster 2 491.204CF+2892.395 | UV-VISo90.449=-0.45TCO-546.715CF+2188.832 

UV-VIS 200-440=0 .781TCO- UV-VIS 200-440=-0 .086TCO- 

Cluster 3 673.848CF+2862.215 635.524CF+2002.077 
U V-VIS200-440= 1.417TCO- 

Cluster 4 1006.816CF+2699.712 | UV-VIS200-440=-0.09TCO-586.346CF+1872.791 
U V-VIS200-440= 1 .603TCO- 

Cluster 5 1108.187CF+2618.299 | UV-VISs09-449=0.244TCO-521.511CF+1668.137 
U V-VIS 200-440=2 ۰ 845 TCO- 

Cluster 6 1230.386CF+2200.607 | UV-VIS200-449=0.295TCO-562.612CF+1598.888 
U V-VIS 599.449=3 .975TCO- 

Cluster 7 1295.804CF+1811.524 | UV-VIS309-449=0.182TCO-835.276CF+1686.477 

UV-VISo00-440=5.11TCO- UV-VIS 200-440=-0.34TCO- 

Cluster 8 1684.114CF+1532.337 1078.882CF+1804.393 


© در معادلات فوق واحد ابرناکی صدم. Joly‏ ازون کلی دابسون و واحد تابش تجمعی AS‏ ژول بر مترمربع در روز می‌باشد. 


-۲-٤‏ مقیاس سالانه 


و در شکل ۱ نیز روند تابش فرابنفش. ازون کلی و ابرناکی در هر منطقه. در مقیاس سالانه بر روی نقشه 
نشان داده شده است. بررسی نتایج تغییرات زمانی تابش فرابنفش )200-440 (UV-VIS‏ در مقیاس سالانه نشان داد 
در مناطق VO‏ ۷ و ۸ روند افزایشی معنی‌داری وجود دارد که این چهار خوشه در نوار ساحلی شمال کشور 
واقع شده‌اند (شکل A‏ در pls‏ مناطق (۱ تا ۶) روند معنی‌داری در (UV-VIS200-440)‏ مشاهده نگردید. 


بیشترین ميزان روند در مقیأاس سالانه تابش فرابنفش UV-VIS200.440‏ مربوط به خوشه A‏ است که به ميزان 


YOA‏ جغرافیا و مخاطرات محیطی شمارة اول 


kj/m?‏ ۱/۶ در سال افزایش یافت. کمترین میزان روند در مقیاس سالانه نیز مربوط به خوشه ۵ و به میزان 


kj/m?‏ ۰/۵ در سال افزایش یافت. 


حدول 1- آماره سطح معنی‌داری (آزمون من -کندال) و شیب روند (Sen)‏ میانگین‌های سالیانه تابش. ایرناکی و 


ازون کلی (۱۹۷۹-۲۰۲۰) 


UV-VIS200-440 CF TCO 

p-value Sen's slope p-value Sen's slope p-value Sen’s slope 

Cluster\Test Kj/m’.year Percent/year (DU/year) 
Cluster 1 0.282 -0.206 0.649 0.0 0.010 0.218 
Cluster 2 1.000 0.000 0.288 0.0 0.017 0.202 
Cluster 3 0.633 -0.065 0.186 0.0 0.190 0.142 
Cluster 4 0.401 0.188 0.020 -0.1 0.315 0.111 
Cluster 5 0.006 0.513 < 0.0001 -0.1 0.362 0.106 
Cluster 6 < 0.0001 1.081 0.002 -0.1 0.333 0.108 
Cluster 7 < 0.0001 1.319 0.002 -0.1 0.214 0.104 
Cluster 8 < 0.0001 1.435 0.008 -0.1 0.159 0.097 


بررسی تغییرات زمانی میزان ابرناکی نیز نشان می‌دهد علت روند افزایشی تابش فرابنفش در مناطق V0‏ 
۷ و ۸ احتمالاً مربوط به روند معنی‌دار کاهشی ابرناکی می‌باشد. همچنین در خوشه ٤‏ نیز روند کاهشی در 
سطح ۹۵ درصد معنی‌دار است» که اثر آن بر تغییرات زمانی م4 0۷-۷180 در این خوشه دیده می‌شود. 
به‌طوری که در خوشه ٤‏ با توجه به جدول V‏ تغییرات 200-440 UV-VIS‏ دارای شیب مثبت و افزایشی است اما 
این تغییرات افزایشی در طی دوره موردمطالعه روند معنی‌داری ندارد. روند کاهشی ابرناکی در خوشه‌های ۶ 
۵ ۰ ۷ و ۸ به میزان ۰/۱ درصد به ازای هرسال در کل دوره است. 

نتایج بررسی روند تغییرات میانگین ازون کلی سالانه در طول دوره‌ی آماری EY‏ ساله نشان داد تغییرات 
زمانی میزان ازون کلی در طول دوره آماری در تمام مناطق افزایشی بوده است. اما طبق جدول V‏ تنها در 
منطقه ۱ و ۲ روند معنی‌دار افزایشی مشاهده می‌شود. نتایج پژوهش چن "و همکاران (۲۰۱۹) نیز با نتیجه به 
دست آمده در این پژوهش مبنی بر افزایشی بودن روند ازون کلی در دوره آماری ۰۲۰۰۵۲۰۱۷ سازگاری 
دارد. 

همچنین بررسی سری زمانی ازون کلی نشان داد در ابتدای دوره آماری ۱۹۷۹ تا ۲۰۲۰ تغییرات کاهشی و 
سپس افزایشی شده است. بررسی‌ها نشان داد در ۸ خوشه موردمطالعه از سال ۱۹۷۹ تا ۲۰۰۰ روند کاهشی 
معنی‌داری وجود داشته است (جدول ۷). بیشترین مقدار روند مربوط به خوشه ۸ است که ميزان ازون کلی 
٤‏ دابسون در سال کاهش ah‏ است و از سال ۲۰۰۱ تا ۲۰۲۰ در هر ۸ خوشه روند افزایشی معنی‌داری 


وجود دارد. بیشترین مقدار شیب روند ازون کلی مربوط به خوشه ٤‏ است که به میزان ۵ دابسون در هر 


1 Chen et al 


سال سیزدهم واکاوی چنددهه‌ای تغییرات تابش تجمعی روزانه فرابنفش سطحی در ایران و... Yog‏ 


دانشمند را بر عهده دارد مطابقت دارد. به‌طوری که عنوان شده است که از سال ۰ تا ۰۲۰۱۳ سطح ازون 
کل ys‏ خر وهای wll ge‏ سات سکره agde J‏ 6 درم Al‏ ات ات دل plod‏ 
را ممنوع کرد. این ممنوعیت در سال ۱۹۸۹ اجرا شد. سطح ازون کلی در اواسط دهه ۱۹۹۰ تثبیت شد و در 


دهه ۰ شروع به بهبود کرد (بانرجی و همکاران Meroe‏ 


حدول -V‏ آماره سطح معنی‌داری (آزمون من -کندال) و شیب روند (Sen)‏ سالانه سری زمانی ازون کلی در دو 
دوره زمانی ۱۹۷۹ تا ۲۰۰۰ و ۲۰۰۱ تا ۲۰۲۰ 


۱۹/۹ ۲۷۰۰۱۰ 
Cluster\Test Sen's slope Sen's slope 
p-value DU/year ان‎ DU/year 
Cluster 1 0.001 -0.283 0.001 0.532 
Cluster 2 0.024 -0.268 0.000 0.595 
Cluster 3 < 0.0001 -0.357 0.001 0.586 
Cluster 4 < 0.0001 -0.367 0.002 0.601 
Cluster 5 < 0.0001 -0.391 0.001 0.625 
Cluster 6 < 0.0001 -0.417 0.001 0.608 
Cluster 7 < 0.0001 -0.423 0.003 0.620 
Cluster 8 0.018 -0.446 0.000 0.578 
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شکل *- روند سالانه میانگین تابش فرابنفش ازون IS‏ و ابرناکی در خوشه‌های موردمطالعه 
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۳۳۹ جغرافیا و مخاطرات محیطی شمارة اول 


همچنین با توجه به جدول ۸ ضریب همبستگی بین مقادیر 0۷-۷150040 با ابرناکی و ازون کلی نشان 
می‌دهد که میزان همبستگی داده‌های 0۷-۷1820044 و ابرناکی» به‌حصوص در مناطق ۷ و LA‏ حدود ۰/1۵ 
بیشتر از همبستگی مربوط به 200-440 LUV-VIS‏ ازون کلی است. بیشترین ضریب همبستگی SU pl‏ مربوط به 
خوشه ۸ و کمترین مربوط به خوشه ۲ بود که به ترتیب در مناطقی واقع شده‌اند که بیشترین و کمترین 
ابرناکی را دارد. افزون بر این» معادلات رگرسیون خطی چند متغیره در هر خوشه در مقیاس سالانه در جدول 
A‏ مشاهده می‌شود. از این معادلات می‌توان برای تخمین مقادیر تابش UV-VIS200-440‏ در مقیاس سللانه 
استفاده نمود. در معادلات رگرسیون SU pl‏ بر حسب صدم ازون GIS‏ بر حسب دابسون و تابش فرابنفش بر 
حسب کیلوژول بر مترمربع در روز است. 
جدول ۸- مقایسه ضرایب همبستگی بین تابش فرابنفش 10۷-۷18044 با ازون (TCO) AS‏ و ابرناکی (CF)‏ و 

معادلات رگرسیون خطی در مقیاس سالانه (۱۹۷۹-۲۰۲۰) 


Clusters Correlation Correlation Linear Regression 

UV/ozone UV/cloud 
Cluster 1 -0.069 -0.754 UV-VIS200-4490=0.25TCO-735CF+253 1.57 
Cluster 2 -0.177 -0.672 UV-VIS300-449=0.067TCO-656.924CF+2402.92 
Cluster 3 -0.250 -0.702 ۲۷-۷ 192: ور‎ 5 
Cluster 4 -0.153 -0.688 UV-VIS00.440=0. 1 17TCO-708.424CF+2265.746 
Cluster 5 -0.033 -0.763 UV-VIS 599.449=0.417TCO-758.377CF+21 8 
Cluster 6 0.041 -0.738 UV-VIS300-449=0.762TCO-796.358CF+1955.643 
Cluster 7 0.051 -0.802 UV-VIS300-449=0.7 86TCO-88 1.785CF+1937.141 
Cluster 8 -0.054 -0.806 UV-VIS200-440=0.532TCO-1093.942CF+2049.134 


۵- نتیجه گیری 


بررسی روند در سری زمانی تابش فرابنفش (انتگرال کلی تابش فرابنفش از طول‌موج ۲۰۰ تا 140 نانومتر 
مربوط به دوره آماری (۲۰۲۰-۱۹۷۹) نشان داد میانگین GLU‏ فرابنفش در مقیاس فصلی در فصل تابستان و 
پاییز در نیمه شمالی دارای روند افزایشی معنی‌دار بوده است و متناسب با آن روند ابرناکی در این فصول در 
نیمه شمالی کاهشی بوده است اما معنی‌دار نبوده است. در فصل پاییز در نیمه جنوبی کشور شامل استان‌های 
سیستان بلوچستان. هرمزگان» Gedy‏ فارس. کرمان» خراسان جنوبی. خلیح‌فارس و بخشی از خوزستان 
(حوشه ۱ و ۲) روند کاهشی در میزان تابش فرابنفش oto‏ شد. علت روند کاهشی مناطق یک و دو روند 
افزایشی ازون کلی می‌باشد زیرا طبق جدول ۳ در این دو خوشه روند افزایشی ازون NS‏ شیب بیشتری 
نسبت به سایر مناطق دارد. این نتیجه‌گیری با نتایج سلاگ و همکاران (۲۰۲۰) مبنی بر افزایش غلظت ازون 
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کلی در اعتدالین تطابق دارد. همچنین ابرناکی نیز در این مناطق تغییرات افزایشی داشته است اما روند معنی- 
دار ندارد. 

در مقیاس سالانه نیز در مناطق شمالی کشور تابش فرابنفش 1۷-۷1500440 دارای روند مثبت می‌باشد اما 
در سایر مناطق روند معنی‌دار سالانه مشاهده نشد. این روند مثبت ناشی از روند منفی میزان ابرناکی در این 
مناطق می‌باشد. روند کاهشی ابرناکی در بخش‌های شمالی و شمال غربی کشور (خوشه ۵» ۷۰۱ و aA‏ 
میزان ۰/۱ درصد به ازای هر سال در کل دوره است و بیشترین میزان روند در مقیاس سالانه تابش فرابنفش 
UV-VIS200-140‏ مربوط به خوشه ۸ است که به میزان kj/m?‏ ۱/۶ در سال افزايش یافت. کمترین میزان روند 
معنی‌دار نیز مربوط به خوشه ۵ و به ميزان /ز ۰/۵در سال افزایش یافت. 

نتایج بررسی آمار dle EY‏ داده‌های باز تحلیل ERAS‏ نشان داد مقدار ازون کلی (TCO)‏ از سال ۱۹۷۹ تا 
۰ شاهد روند کاهشی معنی‌داری در هر ۸ خوشه موردمطالعه بوده است که بیشترین مقدار شیب روند 
مربوط به خوشه ۸ است (با ۰/٤‏ دابسون کاهش در سال). پرخلاف حالت قبل» تمام مناطق در سال‌های ۲۰۰۱ 
الی ۲۰۳۰ روند افزایشی معنی‌دار ازون را تجربه نموده‌اند. بیشینه شیب روند افزایشی ازون AS‏ مربوط به 
خوشه ۶ است که به میزان ۰/۸۲۵ دابسون در هر سال افزایش را نشان می‌دهد. 

نتایج این پژوهش نشان داد به‌طورکلی تابش فرابتفش با ابرناکی همبستگی بیشتری نسبت به ازون کلی 
دارد که این نتایج با نتایج مطالعه سمدلی و همکاران (۲۰۱۲) مبنی بر معنی‌دار نبودن همبستگی ازون و تابش 
فرابنفش مطابقت دارد. فانتولکیس و همکاران (۲۰۲۰) نیز نتیجه گرفتند تغییرات درازمدت در UV‏ عمدتاً 
ناشی از تغییرات در غلظت آئروسول‌هاء ابرناکی و آلبیدو سطح است. درحالی که تغییرات در مقدار ازون کلی 
نقش کمتری در تغییرات درازمدت UV‏ دارد. 

ذکر این نکته نیز ضروریست که همبستگی ضعیف ازون کلی با تابش فرابنفش UV-VIS200-440‏ در 
پژوهش حاضر ممکن است به دلیل لحاظ شدن بخشی از طول‌موج‌های مرئی (1۰۰ تا ۶۶۰ نانومتر) در داده- 
های 0۷-۷1520440 حاصل از خروجی BRAS‏ باشد (زیرا تغیبرات ازون کلی تأثیری بر شدت GRE‏ در این 
محدوده طیف مرئی ندارد). 

به‌طورکلی در دوره آماری موردمطالعه به‌جز مناطق خشک و کویری که به‌صورت کلی ابرناکی بسیار کم 
است» در pl.‏ مناطق ابرناکی روند کاهشی داشت و تابش فرابنفش UV-VIS.‏ حتی با وجود روند 
افزایشی ازون IS‏ در دو دهه اخیر از روند افزایشی پیروی نموده است. در این پژوهش تأثیر غلظت ذرات 
آثروسول بر تابش دریافتی در سطح زمین مورد ارزیابی قرار نگرفت. که با توجه به گسترش دو دهه pel‏ 
خشکسالی در کشور و افزایش فراوانی روزهای همراه با گردوغبار, بررسی تأثیرگذاری این پدیده در 
مطالعات بعدی ضرورت دارد. 
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